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Abstract — We consider a spectrum-sharing network where n 
secondary relays are used to increase secondary rate and also 
mitigate interference on the primary by reducing the required 
overall secondary emitted power. We propose a distributed 
relay selection and clustering framework, obtain closed-form 
expressions for the secondary rate, and show that secondary 
rate increases proportionally to logn. Remarkably, this is on 
the same order as the growth rate obtained in the absence 
of a primary system and its imposed constraints. Our results 
show that to maximize the rate, the secondary relays must 
transmit with power proportional to rT 1 (thus the sum of 
relay powers is bounded) and also that the secondary source 
may not operate at its maximum allowable power. The tradeoff 
between the secondary rate and the interference on the primary 
is also characterized, showing that the primary interference can 
be reduced asymptotically to zero as n increases, while still 
maintaining a secondary rate that grows proportionally to log n. 
Finally, to address the rate loss due to half-duplex relaying in 
the secondary, we propose an alternating relay protocol and 
investigate its performance. 

Index Terms — spectrum sharing, cognitive radio, relaying, 
cooperation, relay selection. 



I. Introduction 

Spectrum-sharing [1], [2] allows unlicensed (secondary) 
users to share the spectrum of licensed (primary) users as 
long as the interference caused on the primary is tolerable. 
This problem is often formulated as maximizing the secondary 
rate subject to interference constraints on the primary, or 
as the dual problem of minimizing the interference on the 
primary subject to a fixed rate for the secondary. Thus, 
reducing the interference footprint of the secondary is of 
paramount interest in spectrum sharing. Multihop relaying and 
cooperative communication is known to significantly mitigate 
interference and increase the sum-throughput in many multi- 
user scenarios [3], among others in broadcast channels [4], 
multiple access channels [5] and interference channels [6]. 
This has motivated the use of relays in spectrum sharing 
networks [7]-[13]. 

This paper studies a spectrum sharing network consisting 
of multiple primary nodes and a secondary system with M- 
antenna source and destination, and n half-duplex relays. 
Unlike conventional relay networks [14], [15], the secondary 
relays must not only maximize the secondary rate but also 
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control the interference on the primary, thus new cooperative 
algorithms are called for. To achieve this goal we propose and 
investigate an approach involving amplify-and-forward (AF) 
relaying as well as relay selection. Under the proposed frame- 
work a closed-form expression is derived for the secondary 
rate, showing that it increases as (Mlogn)/2. Furthermore, 
we propose an augmented scheduling algorithm that recovers 
the half-duplex loss and improves the constant factor in the 
throughput growth rate. Finally, we characterize the trade- 
off between the secondary rate and the primary interference, 
showing that the interference on the primary can be reduced 
asymptotically to zero while the secondary rate still grows 
logarithmically with n. Our results suggest that to maximize 
the secondary rate subject to primary interference constraints, 
one must activate a subset of relays that are chosen based on 
their interference profile on the primary, each of the relays 
transmit with power inversely proportional to n, and the 
secondary source must operate at a power level potentially 
below its maximum available power. These outcomes are 
unique to the cognitive relay networks and are distinct from 
the conventional relay networks, e.g., [14]. 

Some of the related work is as follows. Zhang et al. [16] 
studied the secondary power allocation under various power 
and interference constraints. The throughput limits of 
spectrum-sharing broadcast and MAC were analyzed by Li 
and Nosratinia [17]. Recently, relaying in spectrum sharing 
networks has attracted attention. For secondary outage prob- 
ability Zou et al. [12] and Lee et al. [13] proved that the 
relay selection in spectrum-sharing achieved the same diversity 
as conventional relay networks. For decode-and-forward (DF) 
relaying, Mietzner et al. [8] studied power allocation subject 
to a desired secondary rate, and Asghari and Aissa [9] ana- 
lyzed symbol error rate with relay selection. For AF-relaying, 
Li et al. [10] selected a single relay to maximize the secondary 
rate, and Naeem et al. [11] numerically analyzed a greedy relay 
selection algorithm. 

II. System Model 

We consider a spectrum sharing network consists of N p 
primary nodes and a secondary system with an A/-antenna 
source, an M-antenna destination and n single-antenna half- 
duplex relays, as shown in Figure Q] The average interference 
power caused by the secondary on each of the primary nodes 
must be less than 7 [18]. Let H e C Mxn be the channel 
coefficient matrix from the source to the relays, and F £ 
C nxM and G e £nxN p ^ e ^ c h a nnel coefficient matrices 
from the relays to the destination and the primary nodes, 
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Fig. 1. System model. 
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respectively. Denote h p j £ QMxi as ^ c h anne i vector from 
the source to the primary node £, 1 < £ < N p , The source 
has no direct link to the destination, a widely used model [9], 
[11], [14], [19] appropriate for geometries where the relays are 
roughly located in the middle of the source and destination. A 
block-fading model is considered where all entries of H, F, 
G and h p ^ are zero-mean i.i.d. circular symmetric complex 
Gaussian (CAT) with variance <rf, a\, ai and a^ p , respectively. 

The source communicates with the destination via two hops, 
which in general lowers the required transmit power and thus 
reduces the interference on the primary. In the first hop, the 
source sends M independent data streams across M antennas 
with equal power. The relay i receives 



(1) 



where P s is the source transmit power, which must be less than 
a power constraint P s , s E QMxi - ls [[^ Gaussian signals, 
h* £ QixM j s ,-jjg row j Q f jj^ namely the channel vector 
between the relay i and the source, and m is additive noise 
with distribution cW(0, 1). 

In the second hop, a subset of the relays is selected to 
transmit to the destination. We define a random variable Tj 
to indicate whether the relay i is selected (eligible): 




T, 



1, the relay i is eligible 
0, otherwise 



(2) 



No cooperation among the relays is allowed due to their 
distributed nature. Each relay rotates and scales rj by 



P 



E[T,](P s a s 2 + l)' 



(3) 



where P r is the average relay power and 6i is the rotation 
angle, which are designed in the sequel. Therefore, the signal 
transmitted by the relay i is T^rt. 

After the relay forwarding, the received signal vector at the 
destination is 



y = 




(4) 



where D = diag(TiCi, • • • ,T n c n ) is the relay processing 
matrix and w is the equivalent additive noise. The equivalent 



channel matrix H has entries 

n 

[HjfTjq ^ Ti Ci fmi h 



IQ J 



(5) 



where f m i and hi q are [F] m j and [H]^, respectively. 

In this paper, we focus on the effect of the number of relays 
on the secondary rate, i.e., the so-called "scaling laws" for 
the relays in a spectrum-sharing system. Thus, we allow n to 
increase while N p remains bounded. Analysis of scaling laws 
has a long and established history in wireless communications. 
Among the many examples we mention a few, e.g., [14], [20], 
[21]. 

We refer to cross channels between secondary transmitters 
and primary receivers as interference links. We assume the 
destination knows F, D and H, and the relays only know the 
instantaneous channel gains to which they directly connect, 
i.e., hi and the column i of F. The interference (thus the 
channels) from the primary to the secondary is not explicitly 
modeled for brevity, because its impact can be absorbed into 
the noise term w. 

The cross-channel CSI requirements in a TDD system can 
be met by the secondary nodes detecting packets emitted from 
the primary nodes. Otherwise, under the spectrum leasing 
model [22], the primary nodes can be expected to actively 
promote spectrum reuse by transmitting pilots that can be used 
for cross-channel gain estimation. The latter model applies to 
both TDD and FDD. Regarding the precision of cross-channel 
CSI, only the magnitude of the channel gains are needed, and 
the system can be made robust to imperfections in the cross- 
channel CSI to the relays, as shown in subsequent discussions 
(see Remark [TJ. 

III. Spectrum-Sharing with relay selection and 

CLUSTERING 

Relays that have weak interference links but strong sec- 
ondary links are useful for spectrum sharing, while relays 
that produce a strong interference on the primary may do 
more harm than good. Therefore we use relay selection. In 
spectrum sharing, relay selection and allocation of transmit 
powers are coupled through the interference constraint, an 
issue that is not encountered in conventional (non-spectrum 
sharing) relaying. To make the problem tractable, we propose 
a two-step approach: first the allowable interference per relay 
is bounded, leading to the creation of an eligible relay set. 
Then the secondary rate is maximized by selecting appropriate 
relays from among the eligible set and coordinating their 
transmissions in a manner shown in the sequel. 

A. Eligible Relay Selection 

The interference on the primary nodes is controlled by 
activating only the relays with weak interference links. We 
design the relay selection in a distributed manner that does not 
require CSI exchange among the relays. A relay is eligible if 
and only if all of its own interference link gains are less than 
a pre-designed threshold a. So from (ffj 



T t = 



1, 
0, 



\ga\ 2 < a for 
otherwise 



1. 



(6) 
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where \gu\ 2 is the channel gain between the relay i and the 
primary node £. Note that {|<7ft| 2 }£,i are i.i.d. exponentials 
with mean cr 2 , so {Ti\i are i.i.d. Bernoulli random variables 
with success probability 

p= (l-e- a / a p) N ". (7) 

Since each relay determines eligibility based on its own 
interference links, the eligible relay selection is independent 
across the relays. The average interference from the secondary 
system to the primary node I is 

n n 

^ = 2 E [(E^)(E^D] + ^E[|h p ,,| 2 ] (8) 

= !X>[^| 2 |^1]+^ W 

where the factor | is due to the fact that the relays and the 
source only transmit during half of the time. The second equal- 
ity holds since the design of 0$ is independent of interference 
links, as shown soon. Since Ti = 1 implies \gu\ 2 < a W, we 
have 

E[M 2 |T, ; = 1] < / ^—dx = al-e- a '"v{ a + al) 
Jo a p 

(10) 

= /(«)• (11) 

Combining (O and ( fTTT i. we have W, < 7 if a and P r 
satisfy 

nP r f(a) <max(7 r ,0), (12) 

where 7 r = 2 r y—a^ p P s . As long as ( fT2l holds, the interference 
on all the primary nodes is ensured to be less than 7, 
although the relays are selected distributedly. In our two-hop 
communication the source power P s is chosen so that j r > 0, 
and otherwise the secondary rate is zero. 

Remark 1: We briefly discuss CSI uncertainty in the CSI of 
relay cross-channel gains. Denote the (relay) estimated cross 
channel gain as \gu\ 2 - For simplicity, consider \gu\ 2 has the 
same exponential distribution as the true channel gain \gu\ 2 . 
Assume uncertainty can be modeled as an interval, e.g., that 
the true cross-channel gain is in the interval [0, (1 + e)|gtt| 2 ] 
for some known and fixed e. In this case, if a and P r satisfy 

nP r f(a + e) < max(7 r , 0), 

the interference constraints on the primary will still be ensured. 
Since /(•) is an increasing and bounded function, the impact 
of uncertainty e is to reduce the transmit power at the relays. 

B. Distributed Relay Clustering 

The second part of the proposed method aims to maximize 
the secondary rate. Recall that the source and destination 
have M antennas each; the relays are divided correspondingly 
into M groups {Q m , 1 < m < M}, where each group 
of relays aims to provide a virtual pipe between one of the 
source antennas and the corresponding destination antenna. 
This channel-diagonalization approach is reminiscent of [14] 



but requires more sophisticated analysis because the (eligible) 
relay set is random, as shown in the sequel. 

The relay i £ Q m rotates the received signal by 6i such that 

^ fmihim — I fmi \ \ him \ • (13) 

In this case, all the relays in Q m forward the signal sent by the 
source-antenna m coherently to the destination-antenna m. 

Now, the challenge is to decide the assignment of relays 
to the group Q m , for 1 < m < M. We focus on distributed 
methods so that the coordination among relays is reduced. In 
addition, we decouple the relay clustering from the relay selec- 
tion: the relays decide their groups according to their source- 
relay and relay-destination channels but independent of the 
interference links. Therefore, under this framework, 
and {Ti}" =1 are mutually independent. This decoupling allows 
us to leverage existing relaying methods to enhance the 
secondary rate while bounding the primary interference. It also 
greatly simplifies the analysis. 

We shall consider two clustering schemes: 

1) Fixed Clustering: Here, each of the groups has n/M 
relaysQ Subject to this condition, the relays are assigned to 
the groups in a pre-defined manner. Without loss of generality, 
we assume: 

f (m — l)n mn 1 

G m = [*: [ — w ^ + l< l <—,l<m<Mj. 

2) Gain Clustering: In this clustering we have 

G m = ji : \fkm\ > \h iq \, q ^ m, 1 < q < Afj. 

In other words, the groups are assigned based on the relays' 
channel gain to source antennas. A relay (distributedly) de- 
cides to join in the group m if its gain to the m-th source 
antenna is the stronger than any other channel gains. The 
group assignment of relays is independent from each other 
and is also independent of relay eligibility. Note that Q m is no 
longer fixed but depends on the source-relay channels. Because 
all channels are i.i.d., a relay has equal probability of choosing 
any of the groups. Therefore \Q m \ (the cardinality of Q m ) is 
binomially distributed with parameters (n, -p.). 

IV. Secondary Rate in Spectrum-sharing with 
Relays 

We first derive a general closed-form expression for the sec- 
ondary rate under the proposed framework, and then evaluate 
the achievable rate for specific methods. 

A. Calculation of Secondary Rate 

From ©, conditioned on F, D and H, w is a Gaussian 
vector with autocorrelation 

W = I + FDD t F t . (14) 

The secondary rate in the presence of n relays is denoted with 
R n and is given by: 

i?„-ilogdet(^I+^HHtW- 1 ^ (15) 

'We assume the number of relays n is so that n/M is an integer, however, 
this restriction is not essential and can be relaxed [14]. 
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where i is due to the half-duplex relay constraint. 

Now, we find R n for large n. First, from (|5) and (TTTt , the 
entry of H is 



[H] mg 



Br, 



q = m 
q 7^ to ' 



(16) 



where 



p(o*P a + 1) 



^ \ ^ 1 1i I /mi | | ^tm 1 1 



-O 77-1 777, 



P. 



p(a|P s + 1) 



\ ^ ^ ^ /mi ^im ^ 



j6i 



P 



p(a*P s + 1) ^ 



YT^T 4 / Im ^e^. (17) 



The terms in A mm , B mm and C mg are mutually independent, 
because {T,;}" =1 and are independent from each other. 

So we have the following lemma. 



Lemma 1: If mini< m <M \G> 

A 



w.p.l 



pPr 



1 



PPr 



n\l erf P s 



- J] E[f mi h 



oo as n — > oo, we have 

miH^ml] 0) (18) 



0, (19) 



IMC/ 



PPr 



1 



nV <r 2 P s 



(20) 



Proof: The proof follows from [23, Theorem 2.1] 
and [24, Theorem 1.8.D], and is omitted here. ■ 



From Lemma Q] given \Q, 
H 



w.p.l 



oo Vm, we have: 



diag(ai, • • • ,a M ) 0, 



(21) 



where 



PPr 



n f<T?P, 



— ^ Etl/nriHfttol]. (22) 



ieg„ 



The above analysis indicates that H converges to a diagonal 
matrix for large n (with probability 1). We now show that W 
is also diagonalized as n increases. From (fl4l . we have 



71 T P f f* 

[wl = V + 5 



(23) 



where <L 



1 if m = q and 5 mq = if m ^ q. One 



can verify Kolmogorov conditions [24, Theorem 1.8.D], and 
therefore obtain 



[W] 



mq 



1 f ^ P r E[f mi f* q ] 



// \ — P S CT 2 + 1 



, ^ 0, (24) 



where 



Mfmi\% m = q 



0, 



m^q 



Therefore, we have 
W 

n 



diag(6i, • • • ,b M ) ^ 0, 



(25) 



(26) 



where 



b m = 



P-£LiE[|/™| 2 ] , i 



(27) 



n(P s cr 2 + l) n 

From (ED and d26l i. for large n, the end-to-end channel be- 
tween the source and the destination is approximately decou- 
pled into M parallel channels under the proposed framework, 
where the channel coefficient to is a m and the received noise 
has variance b m . The capacity of this parallel channel is 



M 



^I>g(l 



Mbr, 



(28) 



Therefore, it is reasonable to expect that R n ss R for large 
n. After some calculation (omitted for brevity), we obtain the 
following result. 

Theorem 1: Consider a secondary system with an M- 
antenna source, an A/-antenna destination, and n single- 
antenna relays, in the presence of N primary nodes each 
tolerating interference no more than 7. The secondary rate 
satisfies 

R n -R 0, n -> 00, (29) 



under the proposed relay selection and clustering framework. 

B. Achievable Rate under Specific Clustering Schemes 

We apply TheoremQ]to fixed clustering and gain clustering. 

1) Fixed Clustering: In this scheme, \Q m \ = (so 
Lemma[T]is applicable), and \f m .i \ and \hi m \ are i.i.d. Rayleigh 
random variables with mean and as )f* , respectively. 

Therefore, from (ED, a m = 2L Ms/^P^r> for 1 < to < M. 



7Vt7 s (Trf 

AM Y crJP s +l 

Under this clustering, |/ m i| 2 is i.i.d. exponential with mean a\, 
7n 



and we have b m 
am and b m into 



M 



+ i, for 1 < to < M. Substituting 
, R becomes 

np-K 2 a 2 s a\P r P s 



2 g \ ^ \QM^u 2 d P r + n-\a 2 s P s + I)) 



(30) 



From Theorem Q] under fixed clustering, we have: R n — 

R(f) ^ 0. 

2) Gain Clustering: Since \Q m \ is binomially distributed 

with parameters (n,-h), we have \Q m \/n l/M, and 

Lemma [T] is again applicable. Due to the independence of 
, from (l22l . we have 



I / mi | and |/i im 



1 



pP r 



— J2 E[|/™|]E[|/i im |]. (31) 



where E[|/ roi |] = ^ (i.i.d. Rayleigh) and E[|^ m |] = 
maxi< m <M \h m i\, which is the maximum of M i.i.d. 
Rayleigh random variables. We have 



/i/i = E[ max \h im \] 

L Km<M u 



t. — e 



Af-l 



'■' M , <l.r (32) 



V (_!)Af-m-l^-l^ cr s Mr(|) 



m=0 



(M - to) 3 / 2 



(33) 
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Note that /x/, = aa ~^ for M — 1 (no selection is needed), 
which is identical to the fixed clustering. Based on ( f3Tb and 
\g m \/n l/M, we have a m - ^y^'^' 0. 

Under this clustering, b m remains the same as the fixed 
clustering case, since \f m i\ is still i.i.d. Rayleigh for i g Q m , 
Vto. Substituting a m and b m into d28t . we have 



1 



4M3( CT 2 Pj , +n -l ((T 2p s + 1)) 



(34) 



then: JL 



V. Optimal Power Strategy for Spectrum-sharing 
with relays 

In general, one may envision two competing philosophies 
for relay selection: (1) Allow only relays that have extremely 
weak interference links to the primary. Only very few relays 
will qualify but each of them can transmit at high power. (2) 
Allow a large number of relays to be activated. In this case 
the relay powers must be lowered because not all interference 
links are as "good" as the previous case. 

The key question is: which approach is better? Should we 
use a few select relays with excellent interference profiles, 
or more relays operating at lower power? In this section, 
we optimize the threshold a, the relay power P r and the 
source power P s , while bounding the primary interference. 
The results of this section show that in general the balance 
tips in favor of having more eligible relays operating at low 
power. 

A. Optimal Design of a and P r 

Consider a fixed P s . Since a and P r depend on each other 
via ( fT2l . given a the maximum P r is 



nf(a) 



(35) 



Substituting d35} and © into d30]l and (O shows that 
and i?W attain their maxima (as a function of a) at a = «o 
where: 



a a ~ arg max 



(a 2 s P s + l)/(a) 



(36) 



A closed-form solution for a Q is unavailable but numerical 
solution can be easily obtained. Figure [2] shows the optimal 
design of a based on d36l ). For both fixed clustering and 
gain clustering, according to ((36), a Q = 1.7 maximizes the 
secondary rate. 

Now, we characterize the asymptotic behavior of a a , equiv- 
alently the optimal P r . Because d36t is independent of n, the 
optimal threshold a is not a function of n. So from d35l > the 
optimal average transmit powei@ is P r = 0(n _1 ), i.e., there 
exist real constants di,d 2 > so that dxrT 1 < P r < cfen -1 . 
This implies that the secondary system should on average 
allow many relays to operate at low power. One may intuitively 
interpret this result as follows. To comply with the primary 
interference constraints, the sum power of relays must be 
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Fig. 2. Optimal value of selection threshold a under P s = 5,n = 100 



bounded, and by spreading the total power among more relays 
better beamforming gain is achieved via coherent transmission. 

Now, we study the scaling of the secondary rate. Consider 
examples with fixed clustering (Eq. d30i l) and gain clustering 
(Eq. d34l). If P r — for fixed a (not necessarily 

optimal), we have 

R (f) = M-logn+C u R.M = M \ ogn + d + C 2 , (37) 



where 



M 



log 



pir 2 a 2 s a 2 dlr P s 



mAP(a 2 P s + l)(a 2 plr + f(a)) 



andC 2 =log%.One can view Ci as multi-antenna diversity 
gain by selecting over source-relay channels. From dJTi . the 
secondary rate increases as (M log n)/2, which is summarized 
in the next theorem. 

Theorem 2: Consider a secondary system with an M- 
antenna source, an A/-antenna destination, and n single- 
antenna relays, in the presence of N primary nodes each 
tolerating interference no more than 7. For P r — 0(n _1 ) 
and fixed a, the secondary rate satisfies 



w.p.l 



flog 



1. 



(38) 



under the proposed framework with both fixed clustering and 
gain clustering. 

Theorem |2] holds for a broad class of clustering schemes, 
as long as the corresponding a m and b m are bounded but 
non-zero, i.e., the secondary end-to-end equivalent channel is 
diagonalized with probability 1 as n grows. 

Remark 2: It is possible to extend our results to the case 
of peak interference constraint 7. The secondary source will 
manage its instantaneous interference to be smaller than 7 S on 
all primary nodes by adjusting its transmit power according to 
the largest cross-channel gain to the primaries. Then, the sum 
interference from all the relays must be smaller than j r — 
7 — 7 S . Let P r — £/n where £ is a positive constant. The 

2 lt can be shown that Theorem [JJ still holds if P r scales as 0(n — 1 ). 
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instantaneous interference from all the relays to the primary 
node £ is j£ = £53ILi Ti\gn\ 2 /n. This implies that 

|2l W -P-} 



E[T 2 \gn 



1 



for an arbitrary i € {1, . . . , n}, where we have used the fact 
that Tj and <7& for all i have identical distributions. Therefore, 
£ = 7 r (E[Ti | gft | 2 ]) ensures the instantaneous interference 
on all the primary nodes to be smaller than 7 with probability 
1. 

B. Optimal Source Power 

Due to the primary interference constraints, for any chosen 
a the higher the source power P s , the lower the relay power 
P r , and vice versa. From (130b and (l34l . the rate-maximizing 

P s is 



PI = are 



max 



(2 7 - 



0<P s <min(P s ,27) (27- 



!)/(«) 
(39) 

The unique solution of the above optimization problem is 



P s * = min(P ,P s 



(40) 



where Pn. 



P„ = 



a 2 f(a), otherwise: 



27^ 



9 2 



<r 2 sf(a) 



^(2 1 a 2 f(a) + f 2 (a))(a 2 p + 2 1 a 2 ) 



(41) 



Figure [3] demonstrates the optimal source power as a func- 
tion of three channel parameters a 2 p , a 2 and a 2 d . Three curves 
are shown, in each case one parameter varies while the other 
two are held constant (at unity). In this Figure P s = 10, 
7 = 5 and f(a) — 0.8. As the source-primary channels 
become stronger, the source needs to reduce power; otherwise, 
the relay power must decrease to comply with the primary 
interference constraints, which curbs the rate achieved by the 
second hop. If the source-relay channels become stronger, the 
relay-destination links is the bottleneck and the relays need to 
transmit at higher power, thus once again the source needs to 
reduce power. In contrast, when the relay-destination channels 
become better, the source-relay channels are the bottleneck so 
the source needs to increase power. 

C. Asymptotic Reduction of Interference on Primary 

Multiple relays produce opportunities not only to enhance 
the secondary rate but also to reduce the interference on the 
primary. Suppose the interference on the primary nodes to be 
bounded as 7 = 0(n~' 5 ), which goes to zero as n — > 00. 
From (fT2l . it is sufficient to comply with this constraint if 
P r decreases as <d(n~( 1+s ' > ) and P s decreases as <d(n~ s ). 
Substituting P r and P s into the expression of given 
by (f30b and following some order calculation (the analysis 
of is the same thus omitted), we have 



RW = 



M^Mlogn + 0(l) J<± 
0(1) 8>i 



(42) 
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Fig. 3. Optimal source power with 7 = 5, /(») = 0.8 



The above equation characterizes the trade-off between the 
secondary rate and the interference on the primary: the faster 
of the interference reduction, the slower of the rate growth. It 
also shows that the interference on the primary nodes may be 
mitigated (to zero asymptotically), while the secondary rate 
maintains to increase as 9(logn). 

Remark 3: In the above, the allowable interference 7 is 
made to decline as Q(n~ s ), which leads the growth rate to 
decrease linearly in 5. If 7 is reduced more slowly, e.g., 
decreasing as 9(j^^), the secondary rate can increase at a 
rate of 4£ log n. If we try to mitigate the primary interference 
faster than 0(l/ v / n), the secondary rate will not increase 
logarithmically with n. 

VI. Spectrum-sharing with Alternating Relay 
Protocol 

In this section we consider issues raised by the relay half- 
duplex constraint, i.e., limitations that arise because relays 
cannot listen to the source at the same time as they are 
transmitting. When a subset of relays are activated for relaying 
the previously received information, the inactive relays are 
able to listen and receive information from the source, thus 
in principle the source can transmit continually and the half- 
duplex loss can be mitigated. This is the basic idea of spectrum 
sharing with Alternating Relay Protocol, which is the subject 
of this section. 

The protocol consists of L transmission frames, as shown in 
Figure|4] It is assumed the channel coefficient remains constant 
during each frame, but varies independently from frame to 
frame. The source transmits during frames 1 through L — 1, 
and remains silent during frame L. Since the source transmits 
L — 1 data segments during L time intervals, the rate loss 
induced by the half-duplex relaying is a factor of The 
relays are partitioned into two groups Q\ = {1 < i < ^} and 
Q2 = + 1 < i < n}. During even-numbered transmission 
frames a subset of the relays in Qi transmit to the destination, 
while the relays in Q 2 listen to the source. During odd- 
numbered transmission frames, a subset of the relays in Q 2 
transmit, while the relays in Q\ listen. As shown later, each of 
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Fig. 4. Transmission schedule in the alternating relay protocol (ARP) 



the two relay groups asymptotically achieves a rate that grows 
as 4p -^j-i log n, thus the overall system has a rate that grows 
proportionally to M=^logn. Therefore a good part of the 
half-duplex rate loss can be recovered. 

When either group Qi or group Q 2 is in the transmit mode, 
a subset of relays in the corresponding group is selected to 
transmit. A relay is selected (eligible) if its interference links 
satisfy (0, similar to Section IIII-AI The average interference 
power on the primary node £ takes slightly different forms 
depending on whether L is even or odd. When L is even: 



7e 



L 



El 
2 



E 



2L 



E[|3A 



\Tt = 1] 



(43) 



iee 2 



and when L is odd: 



(L - l)P r 
2L 



»eSi 



|T 4 = 1] 



iee 2 



E[|fl«| 2 |Ti = 1] 



(44) 



To comply with the interference constraints on the primary 
nodes, the threshold a and the relay power P r shall satisfy 



where -f L = 
we use the fact that E 



(45) 

7 - 2a 2 sp P s with P s so that j L > 0, and 

\9ti 



nP r f(a) < max(7 L ,0), 

S 

2 \T,, = ll = /(a). Since from 



Section [V] the optimal P r is proportional to n 1 , we let P r 
77/n, and re-write (@5) as 



77/(0:) < max(7 L ,0). 



(46) 



For the Alternating Relay Protocol, relay clustering is ac- 
complished in a manner similar to Section IIII-BI therefore the 

details are omitted. During frame 2k (or 2k+l), let Q^,\ C Q\ 

i 2 \ : 
(or Q m k C Q 2 ) be the set of relays that assists the antenna 

pair m. As long as minm,k,d \Qm \\ °°> tne secondary rate 
will be obtained following the analysis similar to Section [IV] 
Remark 4: At any point in time, it is possible to allow all 
non-transmitting relays to listen to the source, and be eligible 
to transmit in the next frame. This may give some gains, 
however, it also complicates the relay selection by introducing 
dependence between not only interference links but also other 
links such as source-relay and relay-relay links. It may be 



better for a relay even with a small interference on primary 
to remain inactive if it has also a weak channel to destination 
(therefore it cannot help much) but has a strong channel to the 
source (therefore it can listen well for the next round). Thus, 
any gains will come with a loss of elegance and tractability, 
and therefore this approach is not considered in this paper. 

A. A Simple Example: L = 3 

For illustration purposes, we consider L = 3, where Q\ 
(G2) listen to the source during frame 1, and then transmit 
to the destination during frame 2 (frame 3). We assume fixed 
clustering is used with = n/(2M), for 1 < m < M 

and 1 < d < 2. Let (F^) be the channel coefficient matrix 
between the relays in Q4 and the source (the destination), and 
H,. be the channel coefficient matrix between Q\ and Q2 with 
i.i.d. CAf(0,of) entries. 

We now analyze the rate achieved under Alternating Relay 
Protocol. The optimization of the threshold and the source 
power follows in a manner similar to Section [V] and thus is 
omitted here. 

1) Rate Achieved by Q\: After listening to the source at 
frame 1, Q\ relays to the destination at frame 2. At the end of 
frame 2, similar to the received signal at the destination is 



yi 



rFiDiHxSi + FiDini + wi, 



Hi 



(47) 



where si is the signal sent by the source during frame 1, ni 
is the noise forwarded by the group Q\ of relays, Wi is the 
destination noise. For the group Q\ the relay gains are collected 
into the relay processing matrix 



Dx =diag(r lCl ,--- ,T 9 c 9 ), 



(48) 



where Cj is given by ([3]) so that the average relay power 
constraints are satisfied. One can verify that the equivalent 
channel Hi 



Hi w.p.i T 
— — > pil, 

Jn 



(49) 



where p x = zgjp 
lent noise Wi is 



2p +1 ■ The auto-co variance of equiva- 



1 w.p.l 

-Wi — > Ail, 

n 



(50) 



where Ai = 2 (a^pl+i) ^' Therefore, the end-to-end channel 
is diagonalized for large n, and similar to the results in 
Theorem [TJ the rate achieved during frame 2 satisfies: 



flW-Mlog 1 



npn ajajrjP. 



32AP( W 2 d + 2a 2 P s + 2) 



w.p.l 



0. (51) 



2) Rate Achieved by Q 2 : During frame 2, the relays in Q 2 
receive the signal vector: 



r 2 



^H 2 s 2 + H r Di 




HiSi+nj +n 2 , (52) 



where s 2 is the signal sent by the source during frame 2, 
and the second term corresponds to the interference from the 
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transmission of Q\. During frame 3 the relays in Q2 transmit 
to the destination with a processing matrix 



D 2 = diag(T^ +1 C" + i, • • • ,T n c„), 
where, to satisfy the power constraints, for -| + 



Ci 



V 



(53) 

1 < i < n 

(54) 



np(P s a 2 + 77^2/2+1)- 
At the end of frame 3, the received signal at the destination is 
y 2 =F 2 D 2 r 2 + w 2 

= y^F2D 2 H 2 s 2 + -^F 2 D 2 H r D 1 HiSi 



H 2 



F 2 D 2 H r Dini + F 2 D 2 n 2 + w 2 



(55) 



After correctly decoding Si, the destination cancels the 
inter-relay interference^ i.e., the second term in (1551 . After 
eliminating the inter-relay interference, we have an equivalent 
channel: 



y2 = VM H2S2 



w 2 . 



(56) 



Following steps similar to (fJTJ and (I221 l. we have 

(57) 

Note that w 2 is still a 



H 2 w.p.i 
—= — > p 2 I, 



where p 2 = 2Etp ^ _ f __ ?7 _ 
zero-mean Gaussian vector with auto-covariance 

-W 2 = F 2 D 2 H r DiD|HjD^F^ + F 2 D 2 D^F 2 + I. 

n 

In the right hand side of the above equation, we have 

H r DiD|H} = V H r diag(Ti, • • • ,T») Hj 

np{P s aj + 1) 2 

-I. (59) 



(58) 



lU.p. 1 



Therefore, 

where 

A 2 = 



W 2 — > A 2 I, 



9 9 9 



2(P s a 2 + W 2 /2+l) 



2(P^ 2 



1 



(60) 



1. (61) 



Combining d57| i and ( 1581 , the rate achieved by Q2 is 
where 



i?( 2 >-M log (1 + ^1)^0. 



(62) 



The overall rate is given by the following theorem. 

Theorem 3: Consider a secondary system with an M- 
antenna source, an M-antenna destination, and n single- 
antenna relays, in the presence of N primary nodes each 

interference cancellation requires knowledge of H r at the destination, 
however, we note that even without this knowledge it is possible to obtain the 
same scaling of secondary throughput with the number of relays. Intuitively, 
the inter-relay interference is bounded by a constant that is under our control. 
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Fig. 5. Throughput as a function of source power when n = 100, o"2 
a d = a p = 1 



tolerating interference no more than 7. The secondary rate 
satisfies 

R-{R^+RW)/3 W ^0, (63) 

under the Alternating Relaying Protocol with L = 3 and fixed 
clustering. 

From Theorem [3] the growth rate of R is 



R 



w.p.l 



2Mlogn 



(64) 



Remark 5: Theorem [3] can be generalized to an arbitrary 
number of transmission blocks L. For general L we can 
conclude: 

R w.p.i 



(L-l)M 
L 



log n 



1. 



As L increases, the growth rate of R approaches the maximum 
value of Mlogn. 

VII. Numerical Results 



Unless otherwise specified, we use parameters P s = 10, 
M = 2, N = 1, 7 = 5 and <r 2 s = a\ = 1. 

The secondary rates as a function of source transmit power 
are presented by Figure [5] The theoretical rate under various 
P s is calculated according to (f30b and (l34l . Recall that the 
theoretically optimal P s given by (l40l > is obtained by (l30l 
and (l34l . When the source interference links are very weak, 
e.g., a 2 sp = 0.1, maximizing the source power is optimal, 
which is similar to non-spectrum-sharing networks. When the 
source interference links is strong, e.g., a 2 sp = 1,2, unlike 
non-spectrum-sharing networks, the secondary achieves higher 
rate if the source transmit at power lower than the maximum 
value. This is because the source needs to ensure the relays can 
operate with sufficient power, subject to the total interference 
constraints on the primary nodes. 

Figure [6] verifies Theorem Q] under fixed clustering and gain 
clustering. Here, a — 7, a 2 sp = 1, P r is given by d35l l and 
P s = 5, which is almost optimal as shown in Figure [5] The 
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relays 



simulated average rate of R n under two clustering schemes 
are compared to R^' given by (f30b and given by ( l34b . 
respectively, where the results are well matched for modest 
value of n. The secondary rate increases as the interference 
links of relays become weaker (smaller cr^), since the relays 
can transmit at higher power (but the sum relay power is still 
bounded with n). 

Figure [7] illustrates the tradeoff between maximizing sec- 
ondary rate and minimizing interference on the primary. The 
interference power is 7 = h(n)- & with 5 = 0.1 and 0.2, 
respectively. For 5 = 0.2, the interference power decreases 
faster than S = 0.1, while the secondary rate increases more 
slowly. 

The rate of Alternating Relaying Protocol (Theorem |3) is 
shown in Figure [8] We consider a — > 00, where all the 
relays in Q\ and Q2 transmit alternatively. Here, P s = 5 
and P r is determined by d45] l. The simulated rates match the 
theoretic analysis well under modest value of n. As the relay- 
relay channel becomes weaker (smaller of), the inter-relay 
interference is reduced, and thus the secondary rate increases. 



VIII. Conclusion 

This paper studies spectrum sharing networks with dis- 
tributed AF relaying to improve the secondary rate and reduce 
the interference on the primary. In the asymptote of large 
n (number of relays) the optimal power strategy for the 
secondary source and relays was found, achieving a secondary 
rate proportionally to logn. The half-duplex rate loss was 
reduced and the scaling of secondary rate was enhanced by the 
introduction of the Alternating Relay Protocol. The trade-off 
between the secondary rate and the interference on the primary 
was characterized. Finally, our results show that even without 
cross channel information at the secondary, the secondary rate 
can achieve the growth rate log n. 
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